Effects of Geometric Phases in Josephson Junction Arrays 
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We show that the en route vortex velocity dependent part of the Magnus force in a Josephson 
junction array is effectively zero, and predict zero Hall effect in the classical limit. However, geometric 
phases due to the finite superfiuid density at superconductor grains have a profound influence on 
the quantum dynamics of vortices. Subsequently we find rich and complex Hall behaviors analogous 
to the Thouless-Kohmoto-Nightingale-den Nijs effect in the quantum regime. 
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There have been extensive research activities on the 
vortex dynamics in Josephson junction arrays, where 
physical quantities which determine the vortex dynam- 
ics, such as the vortex potential, the effective vortex mass 
and viscosity, are tunable by nanofabrication techniques. 
One area which has started to attract attentions recently 
is the Hall effect in Josephson junction arrays. In a homo- 
geneous superconductor film it is known that the motion 
of a vortex resembles that of an electron in the presence of 
a uniform magnetic field. The counter part of the Lorentz 
force for an electron is the vortex velocity dependent part 
of the Magnus force (hereafter called the transverse force) 
!|J . Naturally, using the analogy for the fractional quan- 
tum Hall effect in semiconductor heteorojunctions, the 
existence of the quantum Hall effect in a 2-d Josephson 
junction array has been argued by various authors [||-[|. 
These proposals not only possibly have opened a new 
practical way to utilize Josephson junction arrays, they 
also have a fundamental physical implication: the real- 
ization of quantum Hall effect in Boson systems, which 
can be used to test out our theoretical understandings. In 
order to put the above attractive proposal on a firm the- 
oretical ground, a serious and thorough scrutiny should 
be conducted. In the present paper we concentrate our 
attention on the role played by the transverse force. We 
have found that extreme cautions should be exercised 
when using the transverse force in Josephson junction 
arrays. As a consequence, the results which we have ob- 
tained are different from those of early work . 

Our main results are the followings. Because of a very 
large energy required, superconducting grains are inac- 
cessible to vortices. The vortex motion in a Josephson 
junction array is confined to the voids (nonsuperconduct- 
ing areas) and tunneling barriers. Therefore the the local 
contact transverse force is effectively zero in the classical 
limit. However, in the quantum dynamics of the vor- 
tices, the transverse 'force' does play an role because of 
the Aharonov-Bohm type Q scattering of the vortices 
by the superconducting grains. The phase of the wave- 
function for a vortex will be influenced by a vector po- 
tential linked to the finite superfiuid density in the su- 



perconductor grains. The vortices should be considered 
as scattered by a periodic array of the Aharonov-Bohm 
type fluxes and we have set up a tight-binding hard-core 
boson model to capture the main feature. A straightfor- 
ward way to find solutions for this boson model is to map 
it onto a fermion problem. Rich quantum Hall behaviors 
are obtained following the work of Thouless, Kohmoto, 
Nightingale, and den Nijs ||, In the following we present 
our analysis leading to above results. Tunneling junctions 
and a square Josephson junction array will be assumed 
in the present paper. Our analysis can be carried over to 
other types of junctions and arrays with necessary mod- 
ifications. 

We first show that the effect of the contact transverse 
force is zero, but nevertheless geometric phases will be as- 
sociated with the motion of vortices in a Josephson junc- 
tion array. We will start from the nonlinear Schrodinger 
Lagrangian for the motion of the superconducting con- 
densate of a type II superconductor in the clean limit 
because it contains all low energy, long wavelength 
dynamics at zero temperature, such as vortex dynam- 
ics and Josephson relations. The nonlinear Schrodinger 
Lagrangian L is 



L = i%i)*i> Vip* ■ Vip-V(tp*,ip) 

2m* 



(1) 



where m* is the Cooper pair mass, twice of the electron 
mass. The Cooper pair wavefunction ip is normalized to 
the half of the superfiuid electron density. The charge 
of a Cooper pair is 2e. The coupling to the electromag- 
netic field can be put in. Since it will not influence our 
demonstration of the absence of the effect of the contact 
transverse force, it will not be written out explicitly. The 
action is S — J dtcPr L(ip*,ip). Variation of the action, 
SS = 0, gives the nonlinear Schrodinger equation for the 
condensate motion. Writing the wavefunction ip as 



ip = y/p exp(i6 



(2) 



then nonlinear potential V is given by V = (p — 
po) 2 /2N(0), where N(0) is the density of states for each 
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spin projection, and po the average Cooper pair density 
0. Now we are ready to consider the motion of a sin- 
gle vortex in a homogeneous superconductor film. The 
corresponding phase 9 is a function of the vortex posi- 
tion r v (t), satisfying V x V#(r,r„) = 2nq S 2 (r — r v )z. 
Here q = ±1 is the sign of vorticity and z the unit vec- 
tor perpendicular to the film. We consider the vortex 
as a point particle, and r v (t) is defined as the center of 
the vortex core. To obtain a Lagrangian L v for the vor- 
tex, we will perform the integration over r in the action 
S = J dtd 2 r L(r, t, r v ). Only the first term in the nonlin- 
ear Schodinger Lagrangian, Eq.(l), is relevant. We begin 
by 8(r — r„) = V r „#(r — r„) • r v . After the spatial integra- 
tion, the first term in the Lagrangian L in Eq.(l) gives 
rise to a term in the vortex Lagrangian L v as qA. v ■ r v 
with 



qA v 



h / d 2 rp V r „6»(r-r. u ) 



(3) 



Comparing with the known Lagrangian for an electron 
in magnetic field L e = |r 2 + e r ■ A, we conclude that 
the vortex is moving in a fictitious magnetic field. The 
'magnetic field' for a vortex is identified as 



B, 



Vr„ X 



-hpo I . 



(4) 



Therefore the transverse force on a moving vortex is given 
by F m =qr v xB v . 

Next we consider a case in which the Cooper pair den- 
sity is smoothly modulated in space over the scale of the 
vortex core size, p = po(r). In such a case by repeating 
above derivation we have 



B„(r w ) = -q hj d 2 rp(r)V r „ x V r „0(r,r„) 



(5) 



Generally, the phase 9(r,r v ) in a inhomogeneous super- 
conductor consists of two parts: rotational and irrota- 
tional. The irrotational part will not contribute to the 
above integral. The rotational part r still satisfies 



V x V0 r (r,r„) = 2nq 6 2 {r-r v ) 



(6) 



which is required by the single valuedness of the wave- 
function and the presence of a vortex. Using the symme- 
try between r and r' in the above equation we obtain 



B„(r„) = -hz J d 2 rp(r)S 2 (r - r„) = -hp (r v )z 



(7) 



We notice that the fictitious magnetic field B^ is propo- 
tional to the Cooper pair density at the vortex position. 
In other words the transverse force on a vortex is a local 
property, just as the Lorentz force on an electron. Here 
the local Cooper pair density po (r) should be understood 
as the average value over a regime much larger than the 



vortex core, but smaller than the length scale for the 
density varying appreciably, the condition of the smooth 
variation. 

Now we consider the extreme inhomogeneous case: a 
Josephson junction array. If the superfluid density at 
some portions of the superconductor film is smoothly re- 
duced to zero, voids are formed. According to Eq.(7) the 
local contact transverse force is zero for a vortex at a void. 
Voids are connected by low superfluid density regions, the 
tunneling barriers. The typical situation in a Josephson 
junction array is that the energy cost for a vortex to be 
at a tunneling barrier is much smaller than the cost at 
a superconductor grain. Therefore vortices are confined 
to move on voids and tunneling barriers, an example of 
the guided vortex motion. Although the transverse force 
is finite on a tunneling barrier, its effect on the classi- 
cal vortex motion is zero, because vortices must move 
along tunneling junctions. To summarize above analysis, 
the local contact transverse force does not play a role in 
vortex dynamics in a Josephson junction array. Subse- 
quently in the classical limit the transverse force can be 
ignored. We note that our above conclusion differs with 



the one in a recent preprint [9| 
constraint has been applied [1( 



where no such energetic 
. The absence of the en 



route transverse force is in agreement with experimental 
observations on the classical vortex dynamics in Joseph- 
son junction arrays at relatively high temperatures and 
with large effective masses such as vortices moving per- 
pendicular to the driving current jllj , no detectable Hall 
effect jl2|], and the straightline balistic motion |0|]. The 
condition for the classical limit will be given at the end 
of the paper. 

We turn to the Hall effect in the quantum regime. Zero 
transverse force itself does not necessarily exclude all pos- 
sible quantum Hall effects. The difference between the 
classical and quantum dynamics lies in the role played by 
the phase, such as the geometric phase associated with 
the vector potential described in Eq.(3), where the vec- 
tor potential can be finite although the 'magnetic field' 
is zero. This is precisely the case of the Aharonov-Bohm 
effect which illustrates that in quantum mechanics poten- 
tials are more fundamental than forces The analogy 
of the Aharonov-Bohm effect has been studied in Ref. 
[ fL4j . In the remaining part of the paper we explore the 
logical consequences of the geometric phases for the quan- 
tum Hall effect in a Josephson junction array. In order 
to reveal the essential physics and to gain the physical 
insight, we will idealize our problem and approximate 
vortices as hard-core bosons. 

To be commensurated with the existence of the vortex 
inaccessible regions and the geometric phases, we con- 
sider the the tight-binding limit of vortex motion. The 
corresponding Hamiltonian may be written as 



2 



H = t 



a,a 

(1,3) 



I 3 



1,3 



(8) 



whore ai is the boson annihilation operator for a vor- 
tex at j-th void, and ( ) stands for the summation 
over nearest neighbors. The phase Aij is defined on 
the links connected the nearest neighbors, and its sum 
around a plaquette is equal to the geometric phase 2-k$q: 
Hpiaquette Aij = 271-00 . A uniform geometric phase will 
be assumed, where the number of 'fluxes' (f>o is the num- 
ber of Cooper pairs on a superconductor grain, which 
may be controlled by a gate voltage. The interaction be- 
tween vortices is described by Vy, which is long range 
and repulsive. We will treat it as a short range repulsive 
interaction for a first approximation. The tunneling ma- 
trix element t is, in terms of the parameters for a Joseph- 
son junction, an order of y/EjEc exp{— 0{\)^/ Ej / Eq}, 
where Ej is the Josephson junction energy and Eq the 
junction charging energy p5j. The energy scale for the 
repulsive interaction is Ej pq] , which is much larger than 
the tunneling matrix element t. Nevertheless with a con- 
siderable amount of energy two vortices can be put on 
one position. This suggests that vortices are hard-core 
bosons. Then we may approximate the vortex problem 
described by Eq.(8) as a hard-core boson problem, an ap- 
proximation has already been implemented in Refs. 

We are ready to discuss the quantum Hall effect of 
the idealized vortex problem. We do this by mapping 
the hard-core boson problem onto a fermion problem by 
attaching odd number of 'fluxes' on each vortex. This is 
a standard procedure. fl6f| The resulting Hamiltonian for 
the fermion problem is 



,i[Aij+Aij] 



(9) 



(U) 



where Cj is the corresponding the fermion annihilation 
operator at the j-th void. The number of statistical fluxes 
4> s at the j-th void satisfys the constrain ip s = — (2m + 

1) < c] Cj >, With J2pl aq uette A lj = which meanS 

that 2m + 1 fluxes have been attached to each vortex. 
We assume that this mapping will give a mean field so- 
lution with an energy gap separated from its excitations. 
Then the statistical fluxes can be adiabatically smeared 
over the lattice and effectively detached from vortices, as 
shown in Ref. jl7). In this case <p s — — (2m + l)n, with n 
is the magnetic flux frustration, the number of vortices 
per plaquette. Then the resulting mean field problem is 
exactly the Harper- Azbel-Wannier-Hofstadter problem, 
where energy gaps do exist. The quantum Hall behav- 
iors of such a problem have been studied in detail by 
Thouless, Kohmoto, Nightingale, and den Nijs ffj. For 
such a system the quantum Hall conductance a^ H is 



with the integer t r 
tion 



the solution of the Diophantine equa- 



r = s r q + t r p 



(11) 



Here the number of fluxes per plaquette <j) — (j)o — 4> s = 
p/q, with p and q coprime, n = r/q, and r, s r , t r integers 
with \t r \ < q/2. Remember that the mapping has gener- 
ated a Chern-Simons term, which has a contribution to 
the Hall conductance as 



'H 



2m + 1 



(12) 



the quantum Hall conductance of the original vortex sys- 
tem is then rajll 



1 



1 

' H 



1 



' H 



(13) 



Converting back into the electric Hall conductance and 
putting back the unit, we find the electric quantum con- 
ductance of the Josephson junction array is 



oh 



4e 2 1 
h a' 



(14) 



H 



H 



(10) 



It is interesting to observe that a boson-to-fermion map- 
ping has been used here to find the incompressible quan- 
tum Hall states, while in the usual fractional quantum 
Hall effect it is the other way around . As well known 
in the previous study of quantum Hall effect [|l8],[l6) for 
a given set of the 'flux' 0o and the frustration n, there 
may exist several values of m, that is, several mappings, 
with their mean-field solutions all corresponding to filled 
bands which are separated from excitations by energy 
gaps. If such a case occurs, detailed calculation is needed 
fo find the m with the largest energy gap, which is the 
most stable one. 

One can check that following symmetries hold for 
the quantum Hall conductance an- the periodicity, 
<7f/(0o,fi) = cr?/(0o + 1,ti) ; the odd symmetry, 
<JH{4'0i n ) = ~ a H (~ 00) n ); the particle-hole symmetry, 
&H{4>0i n ) = ~ a H (<Poi 1 ~ n )- Because is a non- 
monotonic and rapidly varying function of the 'flux' 
(number of Cooper pairs per plaquette) <j>o and the frus- 
tration n, so will be an- Particularly, both positive 
and negative Hall conductance may be easily reached. 
For example, for n — 1/5 and </>o = 1/3, we find that 
(Th = "y^f - with m = 1; and for n — 1/3 and 0o = 1/5, 
(jff = — |^f- with m = — 1. This is in sharp contrast 
with the previous proposal of the quantum Hall effect 
in a Josephson junction array J|-Q]. There are special 
sets of 4>o and n such that 2m n = <f> , that is, in a 
boson-to-boson mapping the statistical flux cancels the 
real flux. In this case the mean-field solution is automat- 
ically within a gap, and the Hall conductance is 



3 



„ 4e 2 
an = 2m — — 



(15) 



This is to the contrast with the fermion-to-fermion map- 
ping case discussed in Ref. [^0| , where there is no energy 
gap which separate the mean-field solution and excita- 
tions. With these specific sets of 4>o and n and in the 
zero limit of their fraction parts one can take the con- 
tinuous limit of the tight-binding model [|l6| to recover 
previous proposed quantum Hall effect based the contact 
transverse force ]2|-|| . 

Two comments concerning previous work are also in 
order. First, the conditions to observe the Aharonov- 
Casher effect in a Josephson junction array have been 
discussed in Ref. |Q. Because of the extensive spread- 
ing of the associated magnetic field, they are unlikely 
to be fulfilled, which applies to the situation considered 
in Ref. Second, in Ref. jl| the quantum Hall effect 
was argued from the dynamics of Cooper pairs (charged 
bosons), a dual picture of the vortices. Our comment 
here is that even in this dual picture the contact Lorentz 
force is not responsible for the quantum Hall effect, but 
the real Aharonov-Bohm phase similar to the vortex pic- 
ture as discussed in the proceeding paragraphs is. 

To summarize, the effect of the vortex velocity de- 
pendent part of the Magnus force is found to be zero 
in the vortex dynamics in a Josephson junction array. 
Instead the geometric phases are important in the vor- 
tex quantum dynamics. We predict that in the classical 
limit there is no Hall effect at all, but in the quantum 
limit rich quantum Hall behaviors should exist: posi- 
tive and negative Hall conductances determined by the 
Thouless-Kohmoto-Nightingale-den Nijs effect and the 
Chern-Simons contribution. The relevant energy scale 
is the tunneling matrix element t. For a Josephson junc- 
tion energy Ej ~ 1 K, and the junction charging energy 
Ec with a value such that Ej/Ec ~ 1, we find t ~ 100 
mK. When the temperature is higher than t, thermal 
fluctuation will destroy the quantum coherence and the 
vortices move classically. The quantum regime is realized 
for temperatures lower than t where the phase coherence 
is preserved. Experimentally the quantum regime should 
be accessible in principle and a recent report has already 
shown the existence of the Hall effect at low temperatures 
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